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PREFACE 


This  report  was  sponsored  by  the  St.  Paul  District,  US  Army  Corps  of 
Engineers,  under  DA  Form  2544,  number  NCS-IA-85-44-ED-GH  dated  28  Jan  85,  and 
was  monitored  by  Mr.  Dennis  Holme,  St.  Paul  District.  The  report  describes 
the  results  of  calibration  simulations  of  the  water  quality  model  CE-QUAL-R1 , 
representing  Lake  Ashtabula  Reservoir. 

The  model  was  calibrated  and  the  report  was  written  by  Dr.  Joseph 
Wlosinski,  Water  Quality  Modeling  Group  (WQMG),  Ecosystem  Research  and 
Simulation  Division  (ERSD),  of  the  Environmental  Laboratory  (EL),  US  Army 
Engineer  Waterways  Experiment  Station  (WES).  Other  members  of  the  WQMG  who 
assisted  with  the  project  were  Ms.  Sandra  Berry,  Ms.  Dorothy  Hamlin, 

Ms.  Dollie  Sue  Hull,  Mr.  Issac  Jefferson,  and  Mr.  Craig  Oldham.  The  draft 
report  was  reviewed  by  Dr.  Marc  Zimmerman,  Dr.  James  Martin,  and  Ms.  Berry  of 
the  WQMG. 

The  study  was  conducted  under  the  direct  supervision  of  Mr.  Mark  Dortch, 
Chief,  WQMG;  and  under  the  general  supervision  of  Mr.  Donald  Robey,  Chief, 
ERSD;  and  Dr.  John  Harrison,  Chief,  EL.  Director  of  WES  during  preparation  of 
this  report  was  COL  Allen  F.  Grura,  USA.  Technical  Director  was  Dr.  Robert  W. 
Whal in. 

This  report  should  be  cited  as  follows: 

Wlosinski,  J.  H.  1986.  "Calibration  of  a  Water  Quality  Model  for 
Lake  Ashtabula,  North  Dakota,"  Miscellaneous  Paper  EL-86-5, 

US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss. 
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PART  I:  INTRODUCTION 


Water  quality  modeling  is  one  method  available  to  managers  to  help  assess 
and  identify  environmental  factors  which  affect  water  quality  conditions. 

With  a  correctly  calibrated  model,  it  should  be  possible  to  evaluate  the 
effects  of  engineering  alternatives  on  water  quality.  Lake  Ashtabula,  North 
Dakota,  has  become  the  subject  of  such  a  water  quality  study  by  the  St.  Paul 
District,  US  Army  Corps  of  Engineers.  Lake  Ashtabula  is  a  multiple  purpose 
reservoir  located  on  the  Sheyenne  River  (Figure  1).  It  is  formed  by  the 
Baldhill  Dam,  which  is  located  about  19  km  north  of  Valley  City,  North  Dakota. 
The  reservoir  is  operated  by  the  US  Army  Corps  of  Engineers  principally  for 
water  supply,  flood  control  and  recreation  use.  Physical  characteristics  of 
the  reservoir  are  included  in  Table  1. 

This  study  grew  out  of  concern  over  present  water  quality  problems  in 
Lake  Ashtabula  and  possible  effects  due  to  planned  projects.  Planned  projects 
include  the  US  Bureau  of  Reclamation's  Garrison  Diversion  and  the  creation  of 
an  outlet  for  Devils  Lake  by  the  US  Corps  of  Engineers.  Devils  Lake,  a  closed 
subbasin  of  the  Sheyenne  River  watershed,  contains  brackish  water.  Present 
water  quality  problems  in  Lake  Ashtabula  include  frequent  algal  blooms  and 
associated  periods  of  oxygen  depletion  which  have  significantly  reduced  the 
fishery  and  recreation  resource  value  of  the  reservoir.  Possible  strategies 
to  improve  water  quality  at  Lake  Ashtabula  include  changing  the  depth  of 
withdrawal,  nutrient  input  concentrations,  surface  elevation,  or  discharge 
operating  plan,  or  building  a  small  upstream  dam  to  act  as  a  sediment  and 


In  1983  the  St.  Paul  District,  using  data  collected  at  Lake  Ashtabula, 
evaluated  CE-THERM-R1,  a  one-dimensional  reservoir  model  that  simulates 
temperature,  total  dissolved  solids  and  suspended  solids  (Holme , Bakke ,  and 
Wlosinski,  1985).  CE-THERM-R1  is  a  sub-model  of  the  larger  water  quality 
model,  CE-QUAL-R1.  In  that  evaluation,  Lake  Ashtabula  was  represented  as  a 
series  of  three  pools  separated  at  two  reservoir  crossings  (Figure  1).  Each 
pool  was  simulated  separately,  with  the  inflows  to  the  second  and  third  pools 
generated  by  the  model  during  simulation  of  respective  upstream  pools. 

Because  the  results  from  that  study  were  satisfactory,  the  St.  Paul  District 
contracted  for  a  similar  evaluation  of  CE-QUAL-R1  using  data  collected  in 
1981.  In  particular,  the  contract  called  for: 

a.  Making  the  thermal  model  (CE-THERM-R1 ) ,  which  was  developed  on  the 
Boeing  computer  system,  operational  on  the  Control  Data  Corporation  Computer 
system  (CDC). 

b.  Calibrating  the  CE-QUAL-R1  water  quality  model  and  evaluating  model 
performance  using  statistical  methods. 

c.  Allowing  for  evaluation  of  three  hypothetical  reservoir  conditions 
including  a  4-foot  pool  rise,  an  upstream  sediment  trap  Impoundment  and  a 
30  to  40  percent  nutrient  load  reduction. 

d.  Making  the  water  quality  model  (CE-QUAL-R1)  operational  on  the  CDC 
system. 

e.  Furnishing  a  report  that  documents  the  results. 


PART  II:  THE  LAKE  ASHTABULA  MODEL 


CE-QUAL-R1 


In  CE-QUAL-R1,  a  reservoir  is  conceptualized  as  a  vertical  series  of 
horizontal  layers  in  which  thermal  energy  and  mass  are  uniformly  distributed 
in  each  layer.  Horizontal  layer  thicknesses  are  variable  and  dependent  on  the 
balance  of  inflowing  and  outflowing  waters.  Variable  layer  thicknesses  permit 
accurate  mass  balancing  and  reduce  numerical  dispersion  during  periods  of 
large  inflow  and  outflow. 

Inflowing  waters  are  distributed  vertically  based  on  density  differences 
so  that  simulations  of  surface  flows,  interflows,  and  underflows  are  possible. 
Water  density  depends  on  temperature  and  concentrations  of  dissolved  and 
suspended  solids.  Outflowing  waters  are  withdrawn  from  layers  based  on 
density  stratification  using  the  selective  withdrawal  algorithms  of  Bohan  and 
Grace  (1973).  Reservoir  outflows  by  port  can  either  be  specified,  or  the  user 
can  invoke  a  subroutine  which  will  choose  port  flows  in  order  to  meet  a 
downstream  temperature  objective. 

The  heat  budget  includes  the  components  of  short  and  long  wave  radiation, 
back  radiation,  reflected  solar  and  atmospheric  radiaton,  evaporative  loss, 
conductive  heat  transfer,  and  gain  or  loss  through  inflow  and  outflow. 

Vertical  transport  of  thermal  energy  and  mass  is  achieved  through  entrainment 
and  turbulent  diffusion.  Entrainment  determines  the  depth  of  the  upper  mixed 
layer  and  the  onset  of  si  rat  if ical ion.  It  is  calculated  from  the  turbulent 


kinetic  energy  influx  generated  by  wind  shear  and  convective  mixing  using  an 
integral  energy  approach  (Johnson  and  Ford,  1981).  Turbulent  diffusion  is  a 


two-way  transport  process  which  incorporates  a  turbulent  or  eddy  diffusion 
coefficient  that  depends  on  wind  speed,  magnitude  of  inflows  and  outflows,  and 
density  stratification. 

Forces  that  directly  affect  constituent  concentrations  are  temperature, 
irradiation,  wind  speed,  inflow  and  outflow  rates,  and  inflowing  and 
outflowing  masses.  The  physical  distribution  of  mass  is  dependent  upon  the 
diffusive  and  convective  processes  described  above  and  on  settling  processes. 
Biological  processes  also  affect  constituent  concentrations.  Photosynthesis, 
dark  respiration,  photorespiration  and  nonpredatory  mortality  influence  algal 
and  macrophyte  mass.  Grazing  by  fish  and  zooplankton  additionally  influence 
algae.  Ingestion,  egeslion,  and  respiration  affect  zooplankton  and  fish 
growth.  Inorganic  compounds  such  as  ammonia-nitogen  (ammonia), 
nitrite-nitrogen  plus  nitrate-nitrogen  (nitrite-nitrate), 

orthophosphate-phosphorus  (phosphorus),  and  silica  are  consumed  and  produced 
as  a  result  of  photosynthesis  and  respiration.  Phosphorus  and  ammonia  are 
adsorbed  to  solids  according  to  a  modified  equation  for  the  Langmuir  isotherm. 
Ammonia  is  also  removed  by  conversion  to  nitrite-nitrate  under  aerobic 
conditions.  Nitrite-nitrate  is  lost  through  denitr if icaton. 

Detritus  depends  on  algal  and  macrophyte  mortality,  ingestion  by  fish  and 
zooplankton,  zooplankton  egestion,  and  settling.  Decomposition  of  detritus 
contributes  mass  to  ammonia,  phosphorus,  and  inorganic  carbon. 

Inflowing  and  initial  concentrations  for  dissolved  organic  matter  (DOM) 
are  divided  into  labile  and  refractory  DOM  compartments.  Photorespiration 
contributes  to  labile  DOM.  Products  from  DOM  decomposition  are  distributed  to 
inorganic  nutrients  and  refractory  DOM. 


A  *'» 


Dissolved  oxygen  concentration  is  of  primary  importance  to  reservoir 
management.  Oxygen  is  produced  by  algal  and  macrophyte  photosynthesis. 

Oxygen  demand  in  CE-QUAL-R1  is  simulated  by  nitrification,  decomposition  of 
organic  compounds  and  sediment,  respiration,  and  oxidation  of  reduced  products 
of  anaerobic  reactions.  Oxygen  may  also  be  gained  or  lost  at  the  air-water 
interface.  Anaerobic  and  aerobic  conditions  resulting  from  changes  in  oxygen 
concentration  drive  many  other  modeled  processes.  If  the  system  becomes 
anaerobic,  decomposition  of  organic  material  stops,  and  phosphorus,  ammonia, 
dissolved  reduced  manganese,  iron  and  sulfide  are  released  from  the  sediments. 
Sediments  release  almost  all  the  anaerobic  compounds  generated  in  CE-QUAL-R1 : 
reduction  and  inflow  account  for  the  remainder.  Reoxygenation  of  the  system 
reverses  these  reactions. 

Total  dissolved  solids  (TDS)  are  simulated  to  obtain  an  approximation  of 
ionic  strength.  Calculations  based  on  the  equilibrium  reactions  of 
bicarbonate,  carbonate,  and  hydroxyl  ions  and  on  ionic  strength  result  in  the 
pH  value  reported  for  each  layer.  This  value  is  then  used  to  calculate  the 
carbon  dioxide  concentration  which  contributes  to  plant  growth  and  diffuses 
across  the  air-water  interface. 

Total  alkalinity  is  simulated  in  CE-QUAL-R1  to  provide  an  indication  of 
the  buffering  capacity  of  the  system.  Alkalinity  is  modeled  as  a  conservative 
substance,  being  only  advecled  and  diffused.  Suspended  solids  influence  both 
the  density  and  light  regimes.  Suspended  solids  are  subjected  to  advection, 
diffusion,  and  settling.  A  more  detailed  description  of  the  final  model  used 
in  this  study  is  available  in  the  revised  CE-QUAL-R1  User's  Manual 
(Environmental  Laboratory,  1986). 


Modeling  Approach 

For  this  application  of  CE-QUAL-R1 ,  Lake  Ashtabula  was  represented  as  a 
series  of  three  pools  separated  at  the  two  reservoir  crossings,  Keyes  Crossing 
and  Ashtabula  Crossing  (Figure  1).  The  reservoir  width  is  reduced  from  549  m 
to  39  m  and  from  533  m  to  50  m  by  the  bridge  embankments  at  the  two  crossings, 
respectively.  Thus,  each  crossing  was  modeled  as  though  it  were  a  dam  with  an 
outlet  structure,  and  the  one-dimensional  assumption  was  considered  applicable 
to  each  pool.  With  this  approach,  violations  of  the  one-dimensional 
assumption  would  not  be  as  severe,  and  the  headwater  pool  could  represent  the 
settling  basin.  The  physical  characteristics  of  the  three  modeled  pools,  on 
the  first  day  of  simulations,  are  given  in  Table  2.  This  modeling  approach 
follows  from  previous  work  of  the  St.  Paul  District  (Holme,  Bakke,  and 
Wlosinski,  1985). 


Model  Evaluation 


Both  graphical  and  statistical  comparisons  were  made  of  predicted  versus 
measured  values.  The  statistic  used  was  the  Reliability  Index  (RI)  of  Leggett 
and  Williams  (1981).  The  RI  is  scale-variant  and  does  not  depend  on  which  one 
of  the  values  being  compared  is  greater  than  the  other.  In  the  case  of 
perfect  prediction,  the  RI  value  would  be  1.0.  If  all  comparisons  differed  bv 
a  factor  or  two,  the  RI  value  would  be  2.0.  An  RI  of  10  signifies  that  values 
are  an  average  of  one  order  of  magnitude  apart.  An  RI  was  calculated  for  each 
variable  for  each  sampling  period  over  all  depths  as  well  as  for  each  variable 
over  depths  and  sampling  periods.  Over  300  comparisons  were  made  for  each 
pool.  Comparisons  were  made  for  temperature,  total  organic  carbon  (TOC), 
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phosphorus,  ammonia,  total  algae,  nitrite-nitrate,  dissolved  oxygen,  pH,  and 
TDS .  For  each  calibration  simulation,  graphical  and  statistical  comparisons 
were  made  for  all  three  pools.  For  each  pool,  only  the  data  from  the  deepest 
station  was  used  for  model  evaluation.  Only  one  value,  for  each  sampling 
station  for  each  day,  was  measured  for  algae.  Because  at  least  two  values  are 
needed  for  calculating  the  RI,  the  one  measured  concentration  was  arbitrarily 
used  at  0.1  and  1.0  meters. 


Data 


Most  of  the  data  required  for  this  evaluation  had  already  been  compiled 
and  were  available  from  the  St.  Paul  District.  These  data  included  initial 
conditions,  driving  variables  (also  termed  boundary  conditions  or  updates), 
and  calibration  data. 

Initial  conditions  and  calibration  data  were  taken  from  data  collected  on 
a  biweekly  basis  at  1  meter  depth  intervals  during  April  through  September 
1981  (US  Geological  Survey,  1982).  Values  for  initial  conditions  for  the 
three  pools  are  included  in  Appendix  A. 

Meteorological  data  were  obtained  from  the  National  Oceanic  and 
Atmospheric  Administration's  station  at  Fargo,  North  Dakota,  110  kilometers  to 
the  east  of  Lake  Ashtabula.  Data  were  averaged  over  24  hours,  the  simulations 
Lime  step.  Discharge,  temperature,  and  constituent  concentrations,  for  inflow 
to  pool  1,  were  obtained  from  daily  records  at  the  Cooperstown  gaging  station 
(IS  Geological  Survey,  1982).  Outflow  concentrations  from  pools  1  and  2  were 
used  as  inflow  values  for  their  respective  downstream  pools.  An  additional 
tributary  representing  Baldhill  Creek  was  included  in  pool  three.  Measured 
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data  for  this  tributary  were  determined  from  the  Dazey  gaging  station  (US 
Geological  Survey,  1982). 

Coefficients  for  power  curves  which  describe  the  physical  characteristics 
of  the  reservoir  were  estimated  from  sedimentation  survey  data  using 
regression  techniques.  Correlation  coefficients  for  the  stage-area 
relationship  were  0.995,  0.995,  and  0.986  for  pools  1,  2,  and  3,  respectively. 
Coefficients  dealing  with  light  penetration  and  mixing  were  initially 
estimated  using  information  supplied  in  the  User's  Manual  (Environmental 
Laboratory,  1986),  and  were  calibrated  using  measured  data  from  pool  1. 


During  calibration  of  CE-THERM-R1  by  the  St.  Paul  District,  initial 
calculations  to  establish  the  pool  water  levels  and  flows  at  the  two  crossing 
sites  were  unsatisfactory  (Holme,  Bakke,  and  Wlosinski,  1985).  The 
discrepancies  resulted  from  not  including  rain  falling  directly  onto  the  lake 
surface.  Because  rainfall  was  not  Included  in  the  original  version  of 
CE-THERM-R1 ,  the  model  was  modified  to  allow  for  rainfall  events.  Rainfall 
values  were  added  to  the  data  sets  as  driving  variables.  Daily  values  of 
rainfall  averaged  over  stations  located  at  Valley  City  and  Cooperstown,  North 
Dakota,  were  used.  Values  were  obtained  from  the  monthly  summaries  of 
climatological  data  supplied  by  the  National  Oceanic  and  Atmospheric 
Administration.  These  same  changes  were  made  for  the  CE-QUAL-R1  simulations. 
The  amount  of  evaporation  during  the  period  of  simulation  was  estimated  by 
using  pan  evaporation  data  from  the  weather  station  at  Carrington,  North 
Dakota,  multiplied  by  a  pan  coefficient  of  0.7  (US  Department  of  Commerce, 
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1968). 
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Initial  estimates  for  the  biological  and  chemical  coefficients  were  i 

i 

obtained  from  the  CE-QUAL-R1  Users  Manual  (Environmental  Laboratory,  1 9R6 )  and 
from  previous  modeling  studies  (Wlosinski  and  Collins,  1985  a,b).  Simulations 

i 

from  all  three  pools  were  used  for  model  evaluation.  Flux  values,  which  are 
rates  of  change  between  variables,  as  well  as  concentration  predictions,  were 
taken  into  account  when  estimating  coefficients  for  the  next  calibration 
simulation.  A  list  of  coefficients  used  in  final  calibration  simulations  is 
provided  in  Appendix  B. 


Computers 

The  initial  calibration  of  CE-THERM-R1  by  the  St.  Paul  District  used  the 
Boeing  Mainstream-EKS  interactive  time-sharing  computer  system.  Because  the 
Corps  of  Engineers  no  longer  maintains  a  contract  with  this  firm,  the  St.  Paul 
District  required  that  both  CE-THERM-R1  and  CE-QUAL-R1  for  Lake  Ashtubula  be 
made  operational  on  the  Control  Data  Corporation  (CDC)  system.  Calibration 
simulations  for  CE-QUAL-R1  were  made  on  a  VAX  11/750,  the  in-house  computer 
for  the  WQMG,  after  which  time  the  models  were  translated  and  tested  on  thp 
CDC  system. 


Command  Files 


Command  files  were  written,  tested,  and  are  provided  to  the  District  as 
an  aid  in  using  both  CE-THERM-R1  and  CE-QUAL-RI.  Command  files  help  to  make 
the  Lake  Ashtubula  models  "user  friendly".  After  reading  the  command  file. 


the  user  need  only  respond  to  the  computer  questions  to  be  able  to  use  the 
models.  Seperate  command  files  have  been  established  for  CE-THERM-R1 


PART  III:  RESULTS  AND  DISCUSSION 


Statistical  results  from  the  final  calibration  of  the  Lake  Ashtubula 
water  quality  model  are  presented  in  Table  3.  Graphs,  comparing  predicted 
values  (solid  line)  to  measured  data  (x's),  are  presented  for  pool  1  (Figure 
2),  pool  2  (Figure  3),  and  pool  3  (Figure  4).  The  average  RI ,  for  the  nine 
variables  for  which  the  measured  data  were  available,  was  1.98,  2.16,  and  2.08 
for  pools  one  through  three,  respectively.  These  values  compare  favorablv 
with  RI  values,  using  the  same  variables,  from  other  studies.  Calibration  of 
CE-QUAL-R1  for  DeGray  Lake  in  Arkansas  yielded  an  average  value  of  2.63 
(Wlosinski  and  Collins,  1985a),  and  for  Eau  Galle  Reservoir  in  Wisconsin,  2.27 
(Wlosinski  and  Collins,  1985b).  The  overall  RI  for  temperature  from  the 
thermal  model  was  1.07,  and  for  total  dissolved  solids  1.10.  The  comparable 
RI  values  from  CE-QUAL-R1  were  1.08  and  1.11,  respectively. 

All  algorithms  and  coefficients  for  the  Lake  Ashtabula  Model  are  the  same 
as  described  in  the  CE-QUAL-R1  User's  Manual  (Environmental  Laboratorv,  1986), 
except  for  the  utilization  of  nitrogen  by  algae  during  photosynthesis.  Two 
forms  of  nitrogen,  ammonia  and  nitrite-nitrate,  are  modeled.  In  the  original 
model,  the  amount  of  nitrogen  utilized  from  either  compartment  was  based  upon 
the  ratio  of  nitrogen  in  that  compartment  compared  to  total  nitrogen.  This 
assumption  continually  led  to  poor  predictions  (RI  values  above  6.0)  for 
nitrite-nitrate.  The  algorithm  was  modified  to  include  a  coefficient 
representing  the  fraction  of  total  nitrogen,  utilized  during  photosysthes is , 
which  was  removed  from  each  of  the  two  nitrogen  compartments.  This 
formulation  gave  better  predictions  and  was  retained  in  the  model.  The  RI  for 
nitrite-nitrate  after  changes  were  made  ranged  from  3.51  to  4.97  for  the  three 
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pools.  The  nitrogen  cycle  may  not  include  other  processes  important  at  Lake 
Ashtabula.  Peterka  (1970)  found  that  80  percent  of  the  standing  crop  of  algae 
was  due  to  Aphanlzomenon,  a  Cyanophyte  (blue  green  algae).  Blue  green  algae 
are  able  to  fix  elemental  nitrogen  for  use  in  photosynthesis,  a  process  which 
is  not  included  in  the  water  quality  model. 

Although  the  overall  RI  for  all  variables  was  considered  satisfactory, 
individual  graphs  of  some  variables  appear  to  show  that  the  model  does  not 
predict  all  of  the  major  dynamics  measured.  For  example,  predictions  of  algae 
in  pool  3  represent  bloom  conditions  in  late  July  and  early  August,  whereas 
very  little  algae  was  measured  at  the  deepest  station  on  these  dates.  In 
part,  this  discrepancy  is  a  result  of  variability  within  each  pool.  An 
example  of  this  variability,  representing  algae  in  pool  3,  is  shown  In 
Figure  5.  The  solid  line  represents  model  predictions,  the  x's  represent 
algal  concentrations  at  the  deepest  station,  and  the  solid  circles  represent 
algal  concentrations  at  other  stations  in  pool  3.  When  data  from  other 
stations  are  considered,  the  algal  predictions  in  July  and  August  appear  more 
reasonable . 

Oxygen  is  usually  considered  the  most  important  variable  when  assessing 
the  overall  water  quality  of  a  reservoir.  Simulation  of  oxygen  was  considered 
satisfactory,  because  the  model  correctly  predicted  slight  stratification  in 
pools  1  and  2,  and  anoxic  conditions  in  pool  3.  It  appears  that  the  lowest 
concentrations  of  oxygen  in  the  hypolimnion  of  pool  3  occurred  immediatelv 


after  an  algal  bloom  with  very  little  temperature  stratification.  High  winds 
between  the  July  28  and  August  12  sampling  periods  were  probably  responsible 
for  the  reoxygenation  of  hypolimnitic  waters.  Worst  case  conditions  for 


oxygen  would  probably  occur  after  an  algal  bloom  accompanied  by  low  wind, 
heavy  cloud  cover,  and  constant  or  rising  water  temperatures. 

Final  predicted  flux  values  also  appeared  reasonable.  Although  no 
measured  data  were  available  for  comparison  with  predicted  values,  a  number  of 
fluxes  were  compared  to  literature  values.  The  sediment  oxygen  demand  (SOD) 
for  pools  1,  2,  and  3  was  0.40,  0.33,  and  0.31  grams  per  square  meter  per  dav, 
respectively.  These  values  are  within  the  range  of  46  literature  values  for 
lakes  and  reservoirs  as  reported  by  Martin,  Effler,  and  Dobi  (1985). 
Phytoplankton  gross  production  was  predicted  to  be  3.9,  4.0,  and  4.7  grams  of 
oxygen  per  square  meter  per  day.  Measured  values  at  Lake  Ashtabula  during 
1966  and  1967  ranged  from  2.3  to  18.2  grams  of  oxygen  ner  square  meter  per  dav 
(Peterka  and  Reid,  1968). 

The  percentage  of  each  negative  and  positive  oxygen  flux  is  presented  in 
Table  4.  Algal  respiration  and  DOM  decay  were  predicted  to  be  the  most 
important  processes  affecting  oxygen  utilization.  In  pool  1  macrophyte 
respiration  was  also  a  large  sink  for  dissolved  oxygen.  Since  the  maioritv  of 
predicted  DOM  was  created  by  algae,  algal  control  appears  to  be  necessarv  in 
order  to  increase  oxygen  concentrations  in  the  hypolimnion.  In  pool  1,  algal 
photosynthesis  was  responsible  for  most  of  the  positive  flux  of  oxygen, 
whereas,  in  pools  2  and  3  most  of  the  oxygen  is  supplied  through  surface 
exchange.  Controlling  algae  in  pool  1  would  not  have  a  negative  effect  on 
oxygen  concentrations,  because  exchange  of  oxygen  at  the  air-water  boundary 
could  replace  the  oxygen  not  created  during  photosynthesis. 


PART  IV:  SUMMARY  AND  CONCLUSIONS 

A  mathematical  model,  representing  Lake  Ashtabula,  was  developed  based  on 
the  one-dimensional  water  quality  model,  CE-QUAL-R1.  The  reservoir  was 
represented  as  three  pools,  separated  at  two  river  crossings.  This  allowed 
simulation  of  both  longitudinal  and  vertical  variation  in  water  quality. 
Graphical  and  statistical  comparisons  were  made  for  over  1200  predicted  versus 
observed  values  which  were  measured  in  1981.  Variables  included  in  the 
evaluation  were  temperature,  TOC,  phosphorus,  ammonia,  total  algae, 
nitrite-nitrate,  dissolved  oxygen,  pH,  and  TDS. 

The  average  RI  for  all  comparisons  was  1.98,  2.16,  and  2.08  for  pools  1, 
2,  and  3,  respectively.  This  compared  favorably  with  RI  values  from  other 
reservoir  studies.  Predicted  flux  values  also  appeared  to  be  reasonable. 

Both  CE-QUAL-R1 ,  and  the  thermal  model  representing  Lake  Ashtabula 
(CE-THERM-R1),  now  reside  on  the  St.  Paul  District  account  of  the  CDC  computer 
system.  Command  files  for  both  models  were  created  to  allow  District 
personnel  to  easily  use  the  models. 
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Table  1 


Description  of  Baldhill  Dam  and  Lake  Ashtabula  Reservoir 

Baldhill  Dam 
Type 
Length 
Crest 
Top  Width 
Maximum  Height 

Freeboard  above  Project  Pool 

Reservoir 


Contributing  Drainage  Area 

4950  km2 

Elevation 

385.9  meters  msl 

Storage 

84,626,000  m3 

Area 

21,449,000  m2 

Average  Depth 

3.9  m 

Length 

43.4  km 

Maximum  Width 

l  .0  km 

Length  of  Shoreline 

125.5  km 

Mean  Annual  Flow 

105,635,000  m3 

(Sheyenne  River  and  Baldhill 

Creek) 

Average  Residence  Time 

292.4  days 

Compacted  earth  fill 
502.9  meters 

Elevation  389.7  meters  msl 
6.0  meters 
18.6  meters 
3.8  meters 


Table  2 


Physical  Characteristics  of  the  Three  Modeled  Pools 


Pool  1 

Pool  2 

Pool  3 

(Upstream) 

(Middle) 

(Downstream) 

Storage 

11,817,000  m3 

15,523,000 

3 

m 

62,174,000  m3 

Area 

5,220,000  m2 

4,200,000  i 

2 

m 

9,720,000  m2 

Length 

6.7  km 

7.6  km 

19.1  km 

Maximum  width 

0.8  km 

1 .0  km 

1.0  km 

Mean  depth 

2.3  ra 

3.7  m 

6.4  m 

Mean  annual  flow 

92,053,00  m3 

92,053,000 

3 

m 

105,635,000  m 

(Sheyenne  River) 


Average  residence 


47  days 


62  days 


215  days 


Table  3 


Reliability  Index  Values  for  the  Three  Pools  of  the 
Lake  Ashtabula  Water  Quality  Model 


Variable 

Pool  1 

Pool  2 

Pool  3 

Temperature 

1.06 

1.09 

1.10 

TOC 

1.40 

1.58 

1.63 

Phosphorus 

2.53 

3.  j.6 

2.22 

Ammon ia 

2.15 

2.05 

1.89 

Algae 

3.80 

3.27 

4.20 

Nitrite-Nitrate 

3.51 

4.97 

4.24 

Dissolved  Oxygen 

1.18 

1.20 

1.31 

pH 

1.04 

1.02 

1.03 

TDS 

1.11 

1.11 

1.10 

Average 


1.98 


2.16 


2.08 


Table  4 


Percentages  of  Negative  and  Positive  Oxygen 
Flux  for  Different  Processes 


Negative  Oxygen  Flux 

algal  respiration 
ammonia  decay 
detritus  decay 
sediment  decay 
zooplankton  respiration 
anaerobic  oxidation 
fish  respiration 
labile  DOM  decay 
outflow 

refractory  DOM  decay 
macrophyte  respiration 


Positive  Oxygen  Flux 
algal  photosynthesis 
inflow 

surface  exchange 


Pool  1 

Pool  2 

Pool 

29.6 

37.9 

44.0 

3.2 

7.6 

8.0 

1.0 

1.5 

.7 

7.2 

6.1 

5.0 

1.3 

2.5 

1.0 

.2 

.  1 

.2 

1.6 

1.2 

1.0 

15.1 

15.6 

14.6 

4.3 

6.7 

2.1 

8.3 

17.5 

19.3 

28.3 

3.3 

4.0 

50.2 

38.1 

39.4 

2.5 

3.8 

1.4 

12.7 

54.9 

55.5 

3 


macrophyte 
phot  osynt  hes is 


34.6 


3.2 


3.7 
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Figure  4  (continued) 


Predicted  concentrations  of  algae  (solid  line)  versus 
from  the  deepest  station  (x's)  and  the  other  stations 
in  Pool  3. 


APPENDIX  A:  INITIAL  VALUES  OF  STATE  VARIABLES 
FOR  LAKE  ASHTABULA  SIMULATIONS 


Pool 

1 

Pool 

2 

Pool 

3 

Var iables 

Units 

Bottom 

Surface 

Bottom 

Surface 

Bottom 

Surface 

Temperature 

deg .  C 

S.5 

8.5 

7.0 

7.0 

7.0 

7.1 

Oxygen 

mg/1 

11.1 

11.4 

9.1 

9.2 

11.1 

11.1 

Algae  (1) 

mg/1 

.09 

.09 

.09 

.09 

.04 

.04 

Algae  (2) 

mg/1 

.4 

.4 

.25 

.25 

.  1 

.  1 

Algae  (3) 

mg/1 

3.9 

3.9 

2.5 

2.5 

1.0 

1.0 

Zooplankton 

mg/1 

.3 

.3 

.3 

.3 

.3 

.3 

Ammon ia-N 

mg/1 

.  140 

.08 

.14 

.14 

.11 

.09 

N02-N  +  N03-N 

mg/1 

.008 

.01 

.001 

.01 

.14 

.12 

P04-P 

mg/1 

.070 

.080 

.05 

.02 

.09 

.07 

Det  r  itus 

mg/1 

1.5 

1.5 

1.5 

1.5 

1.5 

1.5 

Sed iment 

g/m2 

101 . 1 

101.1 

101.1 

101.1 

101.1 

101.1 

Alkal inity 

mg/1 

337. 

337. 

337. 

337. 

337. 

337. 

Total 

U issolved 
Sol  ids 

mg/1 

403. 

401. 

361 . 

388. 

481  . 

477. 

Suspended 

Sol  ids 

mg/1 

20 . 

30. 

20. 

30. 

20. 

30. 

Labile 

Organ  ics 

mg /I 

14.3 

14.3 

13. 

14. 

12. 

16. 

Refractory 
Organ  ics 

mg/ l 

14.3 

14.3 

13. 

14. 

12. 

16. 

PH 

8.7 

8.6 

8.8 

8.8 

8.9 

8.7 

Part  icu  late 
Manganese 

mg/1 

0.0 

0.0 

n.o 

0.0 

0.0 

0.0 

APPENDIX  B:  COEFFICIENT:',  FOR  THE  LAKE  ASHTUBULA  MODELS 


DESCRIPTION 
*DL  =  dimensionless 
PHYSICAL 

Number  of  outlets 
ports 

Number  of  tributaries 
Reservoir  Latitude 

Reservoir  Longitude 

Reservoir  length 

Minimum  layer 
thickness 

Maximum  layer 
thickness 

Area  coefficient 

Area  coefficient 

Width  coefficient 

Width  coefficient 

Port  1 

Elevat  ion 

Vert ical 
d imens ion 

Hor izontal 
d imens ion 


Model 

Values 

Acronym 

Units 

Pool  1 

Pool  2 

Pool  3 

NOUTS 

DL 

2 

2 

2 

NT RIBS 

DL 

1 

1 

2 

XLAT 

decimal 

degrees 

47.2 

47.2 

47.2 

XLON 

decimal 

degrees 

98.1 

98. 1 

cc 

O' 

RLEN 

m 

8000. 

12964. 

32186. 

SDZMIN 

m 

.4 

.4 

.4 

SDZMAX 

m 

1.0 

1 .0 

1  .0 

ACOEF(l) 

DL 

310200. 

2170000. 

2065800 

ACOEF( 2) 

DL 

1.6624 

.3304 

.6484 

WCOEF(l) 

DL 

218.2 

166.5 

261.6 

WCOEF( 2) 

DL 

.533 

.490 

.317 

EL0UT ( 1 ) 

m 

6.0 

7.3 

8.5 

PVDIM( 1 ) 

m 

2.0 

2.0 

1  .5 

PHD IM( 1 ) 

m 

25. 

25. 

36.6 

APPENDIX  B  (continued) 


Model 

Values 

Descr  ipt ion 

Acronym 

Units 

Pool  1 

Pool  2 

Pool  0 

Port  2 

Eievat ion 

EL0UT(2) 

tn 

1.5 

1.5 

2.9 

Vert  leal 
d imens ion 

PVDIM( 2) 

m 

2. 

2.0 

.61 

Hor  izonta 1 
d imens ion 

PHDIM(2) 

m 

10. 

10. 

1  . 1 

Physical  coefficients 

Turbidity  factor 

TURB 

DL 

2.5 

2.5 

2.5 

Wind  coefficient 

AA 

m / (mb-sec) 1 .2E-9 

1.2E-9 

l . 2E-Q 

Wind  coefficient 

BB 

1  /mb 

1.10E-9 

1.10E-9 

1  .  10F.-9 

Sheltering  coefficient 

SHELCF 

DL 

.96 

.96 

.96 

Penet  rat  ive  convect ion 
fraction 

PEFRAC 

DL 

.00 

.00 

.00 

Wind  mixing  coefficient 

CDIFW 

DL 

.000003 

.000001 

.000003 

Advection  mixing 
coefficient 

CDIFF 

DL 

.000006 

.000006 

.000006 

Critical  density 
for  inflow 

CDENS 

g/1 

.5 

.5 

.5 

Extinction  coefficient 

For  water 

EXCO 

1  /m 

.35 

.35 

.35 

For  inorganic  solids 

EXTINS 

1 /m*mg/ 1 

.11 

.  1  1 

.  1  1 

For  organic  solids 

EXT INP 

1 /m*mg/l 

.11 

.11 

.1  1 

Surface  radiation 
f  ract ion 


n> 


.35 


.35 


SURER AC 


35 


APPENDIX  B  (continued) 


Model 

Values 

Descr Ipt ion 

Acronym 

Un  its 

Pool  1 

Pool  2 

Pool  3 

Reaeration  coefficient 
-  oxygen 

DM02 

,  2 
m  /sec 

1.04E-8 

1 .04E-8 

1 .04E-8 

Reaeration  coefficient 
-  C02 

DMC02 

2 

m  /sec 

1 . 63E-9 

1.63E-9 

1.63E-9 

Fraction  of  dead 
algae  to  detritus 

ALDIGO 

DL 

.25 

.25 

.25 

Organics 

Carbon  fraction  of 
dry  weight 

BIOC 

DL 

.46 

.46 

.46 

Nitrogen  fraction 
of  dry  weight 

BION 

DL 

.09 

.09 

.09 

Phosphorus  fraction 
of  dry  weight 

BIOP 

DL 

.0080 

.0080 

.0080 

Phytoplankton 


Gross  production 

TPMAX( 1) 

1/day 

1.0 

1 .0 

1  .0 

rate 

TPMAX( 2) 

1  /day 

.8 

.8 

.8 

TPMAX( 3 ) 

1  /day 

.6 

.  6 

.6 

Settling  rate 

TSETL(l) 

m/day 

.01 

.01 

.01 

TSETL(2) 

m/day 

.15 

.15 

.15 

TSETL( 3) 

m/day 

.20 

.20 

.20 

Half-saturat ion 

Phosphorus 

PS2P04( 1 ) 

mg/1 

.04 

.06 

.05 

PS2P04( 2 ) 

mg/1 

.06 

.06 

.06 

PS2P04( 3) 

mg/1 

.05 

.05 

.05 

B3 


APPENDIX  B  (continued) 


Descr ipt ion 

Model 

Acronym 

Units 

Pool  1 

Values 

Pool  2 

Pool 

Nitrogen 

PS2N( 1 ) 

mg/1 

.30 

.30 

.30 

PS2N(2) 

mg/1 

.08 

.08 

.08 

PS2N( 3) 

mg/1 

.10 

.10 

.  1" 

Carbon 

PS2C02(  1) 

mg/1 

.05 

.05 

.05 

PS2C02( 2) 

mg/1 

.10 

.10 

.10 

PS2C02( 3) 

mg/1 

.10 

.10 

.10 

Light  saturation 

PISAT( 1) 

kcal/m  /hr  80. 

45. 

30. 

level 

PISAT( 2) 

kcal/m^/hr  80. 

45. 

30. 

PISAT( 3) 

kcal/m^/hr  80. 

45. 

30. 

Maximum  excret ion 

TPEXCR( 1) 

1/day 

.03 

.03 

.03 

rate 

TPEXCR( 2) 

1  /day 

.03 

.03 

.03 

TPEXCR(3) 

1  /day 

.03 

.03 

.03 

Maximum  mortality 

TPMORT ( 1 ) 

1  /day 

.03 

.03 

.03 

rate 

TVMORT ( 2 ) 

1/day 

.03 

.03 

.03 

TPMORT( 3) 

1  /day 

.03 

.03 

.03 

Maximum  respiration 

TPRESP(l) 

1/day 

.085 

.085 

.085 

rate 

TPRESP(  2) 

1  /day 

.070 

.070 

.070 

TPRESP(3) 

1  /  day 

.055 

.055 

.055 

Temperature  multipliers 

Low  threshold 

ALG1T1 

deg.C 

lb. 

lb. 

18. 

ALG2T 1 

deg.C 

6. 

b. 

b. 

ALG3T1 

deg.C 

2. 

? , 

2 

Low  opt imum 

ALG1T2 

deg.C 

22. 

22. 

2?. 

ALG2T2 

deg.C 

lb. 

lb. 

lb. 

ALG3T2 

deg.C 

8. 

8. 

8. 

H igh  opt  imum 

ALG1T3 

deg.C 

2b. 

2b. 

2b. 

ALG2T3 

deg.C 

18. 

18. 

18. 

ALG3T3 

deg.C 

10. 

10. 

10. 

High  threshold 

ALGIT4 

deg.C 

32. 

32. 

32. 

ALG2T4 

deg.C 

24. 

24. 

24. 

ALG3T4 

deg.C 

lb. 

16. 

lb. 

B4 
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Model 

Values 

Descript  ion 

Acronym 

Units 

Pool  1 

Pool  2 

Pool  3 

Low  minimum 

ALG1K1  D 

ALG2KI  D 

ALG3K1  D' 

High  minimum 

ALG1K4  D: 

ALG2K4  D: 

ALG3K4  D! 

Nitrogen  fraction 

FRACN  4(1)  D' 

from  ammonia 

for  algae 

FRACN4( 2 )  D1 

growth 

FRACN4( 3)  D1 

Macrophytes 

Maximum 

TPLMAX  1, 

production  rate 

Maximum  dark 

TMRESP  1, 

respiration  rate 

Maximum  excretion 

TMEXCR  l, 

rate 

Maximum  nonpredatory 

TMMORT  1 , 

rate 

Dead  plants  to 

PLDIG01  D] 

dissolved  organics 

Dead  plants  to 

PLDIG02  D1 

det  r itus 

Dead  plants  to 

PLDIG03  01 

sed iment 

Temperature 

TMPMAC  di 

difference  for 

morta l  Ity 

APPENDIX  B  (continued) 


Descr  tpt ion 


Self  shading 
coe  f  f ic ient 

Carbon 

hal f-saturat ion 

Nitrogen 

half-saturat ion 


Model 

Acronym  Units  Pool  1 


EXTINM  (l/m)x(mg/l)  .04 
PLIMC  mg/1  .05 
PLIMN  mg/1  .01 


PLIMN 


Phosphorus 

half-saturat ion 

PLIMP 

mg/1 

.01 

Plant  density 

PLDENS 

g/m3 

10 

Light  saturation 

PLITE 

kcal/m^/hr 

20 

Nutrient  fraction 
from  sediments 

PLERAC 

DL 

.5 

Maximum  depth  for 
plant  growth 

PLNTDEP 

m 

2.' 

Temperature  multipliers 

Low  threshold 

PLTT  1 

deg.C 

2. 

Low  optimum 

PLTT2 

deg.C 

20 

High  optimum 

PLTT  3 

deg.C 

28 

High  threshold 

PLTT  4 

deg.C 

38 

Low  minimum 

PLTK1 

DL 

.1 

High  minimum 

PLTK4 

DL 

.1 

Zooplankton 

Maximum  ingestion  TZMAX 

Maximum  mortality  TZMORT 


I /day  .65 

1  /day  .02 


Values 
Pool  2 


Pool  3 


APPEND [X  B  (continued) 

Model 

Values 

Description  t 

\cronym 

Units 

Pool  1 

Pool  2 

Pool  3 

Ingestion  efficiency 

ZEFFIC 

DL 

.35 

.35 

.35 

Food  preference 

For  algae  1 

PREF(l) 

DL 

.1 

.1 

.1 

For  algae  2 

PREF( 2) 

DL 

.3 

.3 

.3 

For  algae  3 

PREF( 3) 

DL 

.3 

.3 

.3 

For  detritus 

PREF(4) 

DL 

.3 

.3 

.3 

Respiration  rate 

TZRESP 

1  /day 

.10 

.10 

.10 

Feeding  threshold 

ZOOMIN 

mg/1 

.05 

.05 

.05 

Food  Half-saturation 

ZS2P 

mg/1 

.6 

.8 

.6 

Temperature  multipliers 

Low  threshold 

Z00T1 

deg.C 

2. 

2. 

2. 

Low  opt imum 

Z00T2 

deg.C 

12. 

12. 

12. 

High  optimum 

Z00T3 

deg.C 

25. 

25. 

25. 

High  threshold 

Z00T4 

deg.C 

35. 

35. 

35. 

Low  minimum 

ZOOK  1 

DL 

.1 

.1 

.1 

High  minimum 

Z00K4 

DL 

.1 

.1 

.1 

Fish 

Maximum  Ingestion 

TFMAX 

i/day 

.03 

.03 

.03 

Food  Half-saturat ion 

FS2FSH 

mg/l 

2. 

2. 

2. 

B7 
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APPENDIX  B  (continued) 


Description 


noaej. 

Acronym  Units 


Pool  1 


Va] ues 
Pool  2 


Pool  1 


Food  preference 

For  benthos  and 

FPSED 

sed iment 

For  alga  1 

FPALG( 1 ) 

For  alga  2 

FPALG(2) 

For  alga  3 

FPALG( 3) 

For  zooplankton 

FPZOO 

For  detritus 

FPDET 

Ingestion  efficiency 

FEFFIC 

Maximum  mortality 

TFMORT 

Maximum  respiration 

TFRESP 

Temperature  multipliers 

Low  threshold 

FSHIT1 

Low  opt imum 

FSH1T2 

High  optimum 

FSH1T3 

High  threshold 

FSH1T4 

Low  minimum 

FSH1K1 

High  minimum 

FSH1K4 

Decompos it  ion 

Labile  DOM 

TDOMDK 

Ammonia 

TNH3DK 

APPENDIX  B  (continued) 


Model 

Values 

Descript  ion 

Acronym 

Units 

Pool  1 

Pool  2 

Pool  3 

Detr itus 

TDETDK 

1  /day 

.02 

.02 

.02 

Col  if orms 

TCOLDK 

1  /day 

1.4 

1.4 

1.4 

Sed iment 

TSEDDK. 

1/day 

.008 

.008 

.008 

Refractory  DOM 

TRFRDK. 

1  /day 

.005 

.005 

.005 

Labile  to  refractory 
organics 

TDOMRF 

1  /day 

.01 

.01 

.01 

Nitrite-Nitrate 

denitrification 

TN03DK 

1/day 

.01 

.01 

.01 

Temperature  multipliers 

DOM  low  threshold 

D0MT1 

deg.C 

0. 

0. 

0. 

DOM  opt imum 

D0MT2 

deg.C 

28. 

28. 

28. 

DOM  low  minimum 

D0MK1 

DL 

.1 

.1 

.1 

NK3  low  threshold 

NH3T1 

deg.C 

0. 

0. 

0. 

NH3  optimum 

NH3T2 

deg.C 

20. 

20. 

20. 

NH3  low  minumum 

NH3K1 

DL 

.  1 

.  1 

.  1 

N03  low  threshold 

NO  37  l 

deg.C 

0. 

0. 

0. 

N03  optimum 

N03T2 

deg.C 

28. 

28. 

28. 

N03  low  minimum 

N03K1 

DL 

.1 

.1 

.  I 

Detr  itus 

Detrital  settling 

TDSETL 

m/day 

.2 

.2 

.2 

velocity 
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APPENDIX  B  (continued) 


Mode  1 

Values 

Descript  ion 

Acronym 

Un  its 

Pool  1 

Pool  2 

Pool  3 

Temperature  multiplier 

Dower  threshold 

DETT  1 

deg.C 

1. 

1 . 

1  . 

Opt imum  temperature 
for  decomposition 

DETT2 

deg.C 

30. 

30. 

30. 

Low  minimum 

DETK  1 

DL 

.1 

.  1 

.1 

Chem  ical 

Solids  settling 

TSSKTL 

m/day 

.10 

.10 

.10 

P04  adsorpt ion 

ADSRBP 

1  /m 

3.0 

3.0 

3.0 

Nitrogen  adsorption 

ADSRBN 

1  /m 

40. 

40. 

40. 

P04  adsorption 

ADMAXP 

g/g 

.001 

.001 

.001 

Nitrogen  maximum 

ADMAXN 

g/g 

.003 

.003 

.003 

adsorpt ion 


Stoich  iometry 


Oxygen  - 

ammonia  02NH3 

DL 

4.37 

4.57 

4.57 

Oxygen  - 

detritus  02DET 

DL 

1 .4 

1.4 

1.4 

decay 

Oxygen  - 

respiration  02RESP 

DL 

1 . 1 

1.1 

1.1 

Oxygen  - 

photosynthesis  02FAC 

DL 

l  .4 

1.4 

1  .4 

Oxygen  - 

dissolved  02D0M 

DL 

1.4 

1  .4 

1  .4 

organ  ics 

Oxygen  - 

reduced  02MN2 

m. 

.13 

.15 

.15 

manganese 


BIO 


APPENDIX  15  (continuer!) 


/  Model 

I  Description  Acronym  Units  Pool  1 


oxygen  -  reduced  iron 

02FE2 

DE 

.14 

.14 

.14 

Oxygen  -  suit  ide 

02S2 

DE 

2.0 

2.0 

2.0 

Anaerob  ics 

Dissolved  oxygen 
concentrat ion 

OX  YE  I M 

mg/1 

.S 

.S 

.S 

Sed iment 
th ickness 

SEDTHK 

cm 

S. 

S. 

S. 

Manganese  settling 

TMN4ST 

m/day 

.OS 

.OS 

.os 

Manganese  reduct  ion 
rate 

TMN4RK 

1  /day 

.02 

.02 

.02 

Manganese  release  rate 

TMNRKL 

g/m  /day 

.10 

.10 

.10 

Manganese  oxidation 
rate 

TMN20X 

1  /day 

0. 

0. 

0. 

Particulate  iron 
sett l ing 

TFE3ST 

m/day 

.OS 

.os 

.os 

Iron  reduction  rate 

TFE3RK 

1  /day 

.03 

.03 

.03 

Iron  release  rate 

TFEREE 

g/m/day 

.10 

.10 

.10 

Iron  oxidation  rate 

TFE20X 

1  /day 

0. 

0. 

0. 

Sediment  iron  sulfide 
oxidat ion 

T FES AD 

l  /day 

.90 

.90 

.90 

Iron  sulfide  settling 

T FES ST 

m/day 

.SO 

.SO 

.SO 

Iron  sulfide  oxidation 

TFESBD 

1  /day 

•  S 

.S 

•  S 

Sulfate  reduction 

TS04RE 

1  /day 

.OHIO 

.0010 

.0010 

Sulfur  release  rate 

TSREL 

g/m/day 

.000012 

.000012 

.00001 2 

APPENDIX  B  (continued) 


Model 

Description  Acronym  Units  Pool  1 


Sul  tide  oxidation 

TS20XI 

1  /day 

Sulfide  to  iron  sulfide 
reduct  ion 

TS2DK 

1  /day 

Orthophosphate  sediment 

TXP4RE 

g/m/ day 

Ammonia  sediment 
release 

TCNREL 

g/m/day 

APPENDIX  C:  COMMAND,  SOURCE,  OBJECT,  INPUT,  AND  OUTPUT  FILES 
FOR  THE  LAKE  ASHTABULA  MODELS 


Both  the  Lake  Ashtabula  thermal  model(CE-THERM-R] )  and  water  quality  model 
(CE-QUAL-R1)  reside  on  the  St,  Paul  District's  account  of  the  CDC  computer 
system.  After  establishing  communication  the  user  should  type 


BEGIN, ,ASHPR01 

to  invoke  the  command  file  for  the  thermal  model.  The  first  six  characters 
of  the  command  file  are  letters  and  the  last  character  is  a  number.  For 
the  water  quality  model  type 

BEGIN, , ASH PROQ 

where  all  characters  are  letters. 

Both  models  then  ask  the  user  to  respond  to  a  series  of  questions.  A  typical 
terminal  session  is  included  below,  where  the  small  letters  are  printed  by  the 
computer  and  the  capitals  are  typed  by  the  user. 

BEGIN, .ASHPROQ 

enter  1,2,  or  3  for  pool  to  be  simulated?  1 

enter  outflow  (yes  or  no)  for  outflow  prediction?  YES 

enter  statanal  (yes  or  no)  for  statistical  analysis?  YES 


enter  graphics  (yes  or  no)  for  graphics?  YES 
enter  priority  p2,p3,p4,  or  p6?  P4 
graphics  files  being  prepared, 
after  this  run  is  complete,  type: 

beg  in , , ashgrph 

********warnlng  ******** 

*  graphics  can  be  executed  from  a  tektronics  * 

*  4014  or  equivalent  only.  * 

********warning  ******** 

08.32.20.  submit  complete.  jobname  is  ai3zqrh 

A  yes  answer  to  the  second  question  causes  the  model  to  replace  inflow 
concentrations  for  the  next  pool  with  the  predicted  outflow  concentrations 
from  the  pool  that  was  simulated.  The  second  question  is  asked  only  when 
pool  1  or  2  is  simulated.  For  statistical  analysis  or  graphics,  answer 
YES  to  the  proper  question.  Priority  P2  is  low,  or  overnight,  prioritv, 
and  P b  is  used  for  immediate  turnaround.  P4  represents  normal  daily 
turnaround . 

A  number  of  files  are  associated  with  the  Lake  Ashtabula  models.  Besides 
the  thermal  (CE-THERM-Rl )  and  water  quality  (CE-QUAL-R1)  programs,  five  other 
programs  are  included  in  the  modeling  package.  ORAFORM  is  a  data  manipula¬ 
tion  program,  GRAF8  is  the  graphics  program,  STSTIC  contains  the  statistical 


package,  and  QSPASH  and  TSPASH  update  the  inflowing  concent  rat  ions  for  the 
downstream  pool  for  the  water  quality  and  thermal  models,  respectively. 

The  FORTRAN  programs,  or  source  files,  the  object  files,  and  command  files 
to  compile  the  programs  are  listed  in  Table  Cl.  The  only  files  needed  to  run 
either  model  are  the  object  files.  The  other  files  can  be  written  to  tape  to 
lower  the  cost  of  file  storage. 

Input  data  files  are  also  needed,  and  have  been  created,  for  the  two  models. 
Before  the  user  invokes  the  command  file,  he  or  she  should  edit  the  input  file 
Input  files  which  are  stored  on  the  St.  Paul  CDC  computer  account  are  listed 
in  Table  C2. 

A  list  of  the  output  files  created  by  the  model  are  included  in  Table  C3. 

A  number  of  these  files  are  quite  large.  The  user  should  write  the  files  to 
tape  or  delete  them  in  order  to  save  storage  costs. 


TABLE  Cl 


Source,  Object,  and  Command  Files  for  the  Lake  Ashtabula  Models 


Program 

Name 

Source 

File 

Object 

File 

Command 

File 

CE-THERM-R1 

ASHXCQT 

CQT50BJ 

SFCQT 

CE-QUAL-R 1 

ASHXCQR 

CQR50BJ 

SFCQR 

GRAFORM 

GRAFORM 

GRFOBJ 

S KG FORM 

GRAF8 

ASHGRF8 

GRF80BJ 

SFGRF8 

ST STIC 

ST  STIC 

STOBJ 

QSPASH 

QSPASH 

QSPOBJ 

SFQSP 

TSPASH 

TSPASH 

TSPOBJ 

C4 


TABLE  C2 


A  List 

of  Input  Files  for  the 

Lake  Ashtabula 

Mod  els 

File  Name 

Used  in  Program 

Used  for  Pool 

Model 

ASHTP1 

CE-THERM-R1 

1 

thermal 

ASHTP2 

CE-THERM-R1 

2 

thermal 

ASHTP3 

CE-THERM-R1 

3 

thermal 

ASHQP1 

CE-QUAL-R1 

1 

water  quality 

ASHQP2 

CE-QUAL-R1 

2 

water  qual ity 

ASHQP3 

CE-QUAL-R1 

3 

water  quality 

GRP1 

GRAFORM 

1 

thermal 

GRP  2 

GRAFORM 

2 

thermal 

GRP3 

GRAFORM 

3 

thermal 

GQ  PI 

GRAFORM 

1 

water  quality 

GQP2 

GRAFORM 

2 

water  qual ity 

GQP3 

GRAFORM 

3 

water  quality 

G8TP1 

GRAF  8 

1 

thermal 

G8TP2 

GRAFS 

2 

thermal 

G8TP3 

GRAFS 

3 

thermal 

G8QPI 

GRAFS 

1 

water  quality 

G8QP2 

GRAF8 

2 

water  quality 

G8QP3 

GRAFS 

3 

water  quality 

VASHP1 

GRAF8  and  STSTIC 

1 

both 

VASHP2 

GRAFS  and  STSTIC 

2 

both 

VASHP3 

GRAFS  and  STSTIC 

3 

both 

stsvtt 

STSTIC 

1,2,3 

thermal 

STSWQ 

STSTIC 

1,2,3 

water  qualitv 

TABLE  C3 


A  List 

of  the  Output  Flies 

for  the  Lake  Ashtabula 

Model s 

File  Name 

From  Program 

From  Pool 

Model 

PLTWC 

CE-THERM-Rl 

1,2,3 

thermal 

PLT  IN 

CE-THERM-R1 

1,2 

thermal 

PLTOUT 

CE-TUERM-R1 

1,2,3 

thermal 

OFILE 

CE-THERM-R1 

1,2,3 

thermal 

PLT INI 

CE-THERM-R1 

3 

thermal 

PLTIN2 

CE-THERM-R1 

3 

thermal 

PLQWC 

CE-QUAL-R1 

1,2,3 

water  qual  Ity 

PLQIN 

CE-QUAL-R1 

1,2 

water  qual  Ity 

PLQOUT 

CE-QUAL-R1 

1,2,3 

water  quality 

QF 1 LE 

CE-QUAL-R1 

1,2,3 

water  quality 

WCFLOW 

CE-QUAL-R1 

1,2 

water  qual  ity 

PLQIN1 

CE-QUAL-R1 

3 

water  quality 

PLQLN2 

CE-QUAL-R1 

3 

water  qual  ity 

GRAOTP 

GRAFORM 

1  ,2,3 

both 

ST  ST 

ST STIC 

1,2,3 

both 

TAPE  16 


r.RAFR 


1  T  •> 


